Electronic parameters for the hole transfer in DNA duplex oligomers. 
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We report on our calculations of the inner-sphere reorganization energy and the interaction of the 
n orbitals within DNA oligomers. The exponential decrease of the electronic coupling between the 
highest and second highest occupied base orbitals of the intrastrand nucleobases in the (A-T)„ and 
(G-C)„ oligomers have been found with an increase of the sequence number n in the DNA structure. 
We conclude that for realistic estimation of the electronic coupling values between the nucleobases 
within the DNA molecule, a DNA chain containing at least four base pairs is required. We estimate 
the geometry relaxation of the base pairs within the (A-T)„ and (G-C) n oligomers (n = 1 — 6) due 
to their oxidation. The decrease of the inner-sphere reorganization energy with elongation of the 
oligomer structure participating in the oxidation process have been observed. The maximum degree 
of geometry relaxation of the nucleobase structures and correspondingly the higher charge density 
in the oxidized state are found to be located close to the oligomer center. 
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I. INTRODUCTION 

Discovery of charge migration in DNA molecules has 
opened new avenues to investigate various possibilities 
rangingfrom its role in the DNA oxidative damage and 
repair [1] to application of DNA in nanoelectronic de- 
vice developments [2]. In fact, DNA-based molecular 
electronic devices are expected to operate within the pi- 
coseconds range [!, 0| that can exceed the potential of 
the present solid state devices. Quite expectedly, the 
DNA molecule has become a subject of intense research 
activities both theoretically [j, ja HI, H, OH EH, an d 
experimentally @, El, EI M, 

From all these studies of charge migration in the 
DNA molecule reported as yet, it is clear that there 
are two mechanisms for transfer of charge depending on 
the DNA structure and transfer parameters: a superex- 
change charge transfer and the incoherent hopping Q. 
The charge migration leads to the geometry changes in 
the nucleotides and the surrounding environment, which 
significantly contribute to the charge migration process. 
Due to the interaction of the it orbitals of the near- 
est neighbor duplexes and insignificant IP difference be- 
tween them, hole can be distributed over several sites 
in the (A-T)„ and (G-C)„ oligomers. This significantly 
changes the magnitudes of the geometry relaxation of 
the nucleobases - inner-sphere component and environ- 
ment contribution - outer-sphere component. However, 
the investigation of the transfer parameters, such as or- 
bital overlapping 0, El, E3] and activation energy for 
charge migration i.e. the IP and the reorganization en- 
ergy [1, [t| E3, E3, El, Ell HH , nave been performed mostly 
for the nucleobases or/and base pairs. 

The main purpose of our work is to estimate the elec- 
tronic coupling between the two nearest nucleobases, 
their charge distribution and inner-sphere reorganization 
energy, when they are placed within the (A-T)„ and (G- 
G) n oligomer duplexes. All these computations have been 



performed using accurate quantum-chemical methods. 



II. METHOD OF COMPUTATION 

The relatively small reaction free energy in the DNA 
molecule makes the DNA hole transfer mechanism qual- 
itatively different from that in most proteins [22| . The 
electron transfer in the DNA molecule was found to be 
strongly dependent on the details of the donor and ac- 
ceptor energies and deviation of their geometries EH- 

The charge transfer in a DNA molecule occurs due to 
the overlapping between the 7r-electrons of the carbon 
and the nitrogen atoms that forms the ir — ir orbitals be- 
tween the parallel nucleobases. Charge migration in the 
molecular systems with weakly interacting donors and 
acceptors, such as between the base pairs in the DNA 
molecule, is described by the standard high-temperature 
nonadiabatic electron-transfer rate 



k = ^-\H DA \ 2 (FC), 



(1) 



where Hda is the electronic donor-acceptor matrix ele- 
ment, and FC is the Franck-Condon factor. 

The electronic donor-acceptor matrix element Hba is 
defined by the coupling of the orbitals of the donor and 
the acceptor and depends on the structure of the DNA 
molecule. For the (A-T)„ and (G-C) n oligomers the sim- 
ple expressions for the deviation of the electronic coupling 
on the sequence number n have been generated [23[ . Ac- 
cording to these expressions, the value of the electronic 
coupling decreases with elongation of the oligomers [23[ . 
In Sect. Ill A, we simulate the electronic coupling of the 
nucleobases within the (A-T)„ and (G-C)„ oligomers us- 
ing the quantum chemistry methods with the Jaguar 6.5 
program [24[ . According to the Koopmans' theorem, the 
electronic coupling can be estimated as half of the adia- 
batic state splitting between the HOMO and the HOMO- 
1 of the closed shell neutral system, determined in a 
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FIG. 1: Schematic diagram for the Marcus-Hush scheme of 
the inner-sphere reorganization energy: (a) the reorganization 
energy Xi : o for adding one electron to the positive ion and (b) 
the reorganization energy A^a for removing one electron from 
the neutral geometry, where D and A are the hole donor and 
the hole acceptor. 

Hartree-Fock self-consistent field. Therefore, the RHF /6- 
31 + G* have been applied for the electronic coupling cal- 
culations. The 6-31 + G* basis set is appropriate for our 
purposes. Previous investigations indicated that any fur- 
ther extension has little influence on the electron coupling 
@ . The geometries of the separated DNA base pairs have 
been optimized with the RHF/6-3I + G* bases and in the 
following, the optimized geometries of the base pairs have 
been stacked with a twist angle 36° and a distance of 3.38 
A. The stacking of the preliminary optimized geometries 
allows us to consider the same nucleobases to be 'in reso- 
nance' [9( within the structures of the (A-T) n and (G-C)„ 
oligomers. 

The FC factor deals with the influence of the vibronic 
interaction on the charge propagation and can be ex- 
pressed as 

TO -(«*■*,-» «p(-^±£), (2) 

where AG is the free energy of the reaction, and A = 
X s + Xi is the reorganization energy. The interaction of 
the molecule with the solvent environment is included 
in the outer-sphere reorganization energy A s , while the 
relaxation of the acceptor, the donor and the molecular 
bridge geometries are included in the inner-sphere reor- 
ganization energy A^. 

The inner-sphere reorganization energy accounts for 
the low-frequency inter-molecular modes and can be es- 
timated within the quantum chemical approach as (25| 

Xi = Xi t A + Xi t D = (E+(A) ~ E +(A)) + (E( D ) -Sp)), (3) 

where E is the energy of the neutral state in a neutral 
geometry, E* is the energy of the neutral state in an 



ionic geometry, E+ is the energy of the ionic state in an 
ionic geometry, and El. is the energy of the ionic state 
in a neutral geometry. The reorganization energy X^a 
is the energy to remove an electron from the hole accep- 
tor A, while the reorganization energy X^d is the energy 
to add an electron to the hole donor D + . The scheme 
for calculation of the reorganization energy is presented 
in Figure [T] Clearly, the vibronic interactions stabilize 
the geometry of the donor and the acceptor from a non- 
equilibrium state (E^_, A s, E? D s) to the equilibrium state 
(E+(m, E(m). The vertical ionization potential is deter- 
mined as vIP = (J5* — E) and differ from the adiabatic 
IP = (E + — E) by the inner-sphere reorganization energy. 

The inner-sphere reorganization energy has been eval- 
uated within the unrestricted Becke3P86/6-31f + G* ap- 
proximation of the DFT method. The DFT theory was 
found to be reasonable for this purpose based on a com- 
parison of the results of Ref. [2I(. These results show 
that the DFT theory predicts the magnitude of the inner- 
sphere reorganization energy with a minimum error when 
compared to the experimental data [22|, HH, H3- Fur- 
thermore, we have also tested the application of the HF 
method and the DFT theory for the vertical ionization 
potential (vIP) calculations and have found significant 
qualitative and quantitative disagreement of the HF with 
the experimental data [22j], while the Becke3P86 approx- 
imation is appropriate for this purpose. 

III. RESULTS AND DISCUSSION 
A. High occupied orbital distribution 

At first we consider the system of two stacked du- 
plexes. The results for the highest occupied base or- 
bital (HOBO) are presented in Table fl] In the case 
when the pyrimidine/pyrimidine and purine/purine bases 
are stacked in one strands, the HOBOs of the adenine 
and guanine bases have lowest energy in comparison to 
the pyrimidinc/purine configurations. For the (A-T)2 
and (G-C)a oligomers the HOBOs are delocalized over 
the two intrastrand nucleobases, and therefore, it pro- 
duces a significant coupling between the 7r orbitals of the 
stacked pyrimidine/pyrimidine and purine/purine bases. 
For oligomers where the pyrimidine and the purine bases 
are stacked in the same strand (A-T/T-A, G-C/C-G, A- 
T/C-C) or in the mixed structures (A-T/G-C and G- 
C/A-T), for some cases the ir orbitals are delocalized, but 
electronic coupling is weak. For others the 7r orbitals are 
localized mostly on one nucleobase, and we can consider 
the weak intrastrand and interstrand coupling between 
the nucleobases as well. The low interstrand coupling 
for the cases A-T/T-A and G-C/C-G has been observed 
experimentally [28[. 

In the (A-T) n and (G-C)„ DNA oligomers, the 2n tt 
orbitals are delocalized over the neighboring intrastrand 
nucleobases. We have found that the HOMO corre- 
sponds to the central nucleobases in the (A-T) n and (G- 
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TABLE I: The HOBO within the system of two stacked base pairs estimated with RHF/6-31+G7/RHF/6-31+G*, in the case 
when the orbital is localized (1) on a single nucleobase, or delocalized (d). All values are in eV. 
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9.33(1),9.36(1) 


9.71(d) 


9.81(d) 


9.65(1) 



C)„ oligomers structures. Particularly, the HOBOs of 
the adenine and guanine primarily belongs to the nu- 
cleobase in the center of the oligomer, while HOBO-1 
belongs to the nearest neighboring nucleobase to that in 
the oligomer center. Similarly, the HOMO resides pri- 
marily on the central guanine have been observed for the 
5'-(G-C)„ structures [2{|. This effect is related to the 
electrostatic potential distribution over the (A-T)„ and 
(G-C)„ oligomers in the vacuum. We have computed the 
electrostatic potential distribution in the (A-T)„ and (G- 
C) n oligomers with the the APBS program performing 
the nonlinear Poisson-Boltzmann solver [3(| . The RESP 
procedure [3l| has been applied to determine an atomic 
partial charge of the A-T and G-C base pairs. The atomic 
partial charges of the base pairs have been the same for 
each base pair in the oligomer structures. We have found 
that for the (A-T)„ structure, the the electrostatic poten- 
tial on the central adenine and thymine is more negative 
than on the nucleobases of the oligomer sides. For the 
(G-C) ra structure, the electrostatic potential on central 
guanine is more negative than on the sides, while for the 
cytosine the opposite effect takes place. Therefore, the 
stronger localization of the HOMO density on the central 
nucleobase is observed for the (A-T)„ oligomer, where the 
HOMO is mostly delocalized over three neighboring in- 
trastrand adenines, than for the (G-C)„ oligomer, where 
the HOMO is delocalized over four neighboring guanines. 
The HOMO electronic density for the (A-T) 5 DNA se- 
quence are presented in Figure [U It has been found, 
that for the n=3 the population of the HOMO is much 
larger for the central nucleobases than for the sides. For 
the n > 3 the population of the HOMO on the central 
nucleobases decreases due to the HOMO dclocalization. 

The derealization of the orbital electron density over 
the oligomer structure produces a decrease of the HOBO 
energies with elongation of the DNA chain. The depen- 
dence of the orbital energies of the n HOBOs, which are 
7r orbitals of the nucleobases, in the (A-T)„ and (G-C) ra 
oligomers on the sequence number n are presented in Fig- 
ure [3l As we see from the results, the splitting of the 
HOBO and HOBO-1 decreases with elongation of the du- 
plex oligomer structures. We conclude that the electronic 
coupling between the HOBO and HOBO-1 belonging to 
the nearest intrastrand nucleobases decreases as well due 
to the spreading of the electron density of the molecular 
7r orbitals over the larger sequence number. Therefore, 




FIG. 2: The electronic density of the HOMO that resides 
on the central adenine (as indicated) in the (A-T)s DNA se- 
quence (RHF/6-31 + G*//RHF/6-31 + G*). 



the maximum value of the electronic coupling is observed 
for the structures of the two stacked base pairs, where the 
HOBO can spread only over two nucleobases. The elec- 
tronic coupling between the two HOBO and HOBO-1 for 
the different nucleobases within the (A-T) n and (G-C) n 
oligomers are presented in Figure 2] As we mentioned 
above, the ir orbitals are primarily delocalized over 3- 
4 intrastrand nucleobases. Therefore, a fast decrease of 
the electron coupling for n <4 is observed when the 7r 
orbitals have the potential to spread (see Figure \2§ ■ For 
n > 4 the electronic coupling decreases slowly. For n=6, 
the electronic coupling magnitude is less than half of that 
for n=2. According to the performed extrapolation pro- 
cedure in Figure [H for n >8 the coupling is practically 
independent of the sequence number. 

The observed decrease of electronic coupling in the (A- 
T)„ and (G-C) n oligomers is in good agreement with the 
approximation in Ref. [23j |. where for n=4 the electronic 
coupling is half of that for n=2. Moreover, it can ex- 
plain the disagreement of the earlier theoretical results 
[tj, [13] and experiments, where the electronic coupling is 
usually much smaller than 0.1 eV [HI]. Based on these 
results we conclude that the electronic coupling calcu- 
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TABLE II: The inner-sphere reorganization energy and the 
values of vIP for the nucleobases and pairs estimated with 
UB3P86/6-311+G*. All values are in eV. 



1 2 3 4 5 
sequence number (n) ► 



FIG. 3: The n orbitals of the nucleobases within the (A-T) n 
and the (G-C) n oligomers, where n is number of the DNA 
base pair (RHF/6-31+G*// RHF/6-31+G*). 
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FIG. 4: The electronic coupling between the HOBO and 
HOBO-1 for the nucleobases in the (A-T)„ and the (G- 
C)„ duplex oligomers (RHF/6-31+G*// RHF/6-31+G*). The 
dotted lines correspond to the extrapolated results. 



lated for the system of two stacked base pairs leads to 
the overestimation of the transfer integral for the charge 
transfer simulation in the DNA oligomers. It seems to us 
that more accurate value of the electron coupling can be 
determined from the DNA chain with length of at least 
four base pairs. 



B. Inner-sphere reorganization energy 

The orbital overlapping and the distribution of the HO- 
MOs over the oligomer structure directly influence the 
charge distribution in the DNA molecule and the inner- 
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sphere reorganization energy, respectively. 

We have calculated the inner-sphere reorganization en- 
ergy and the vIP for the separated nucleobases and their 
pairs. The results are compared with the experimental 
data in Tablc|TTl From Table HJ it is clear that our calcula- 
tion and the experimental results for the nucleobases are 
in good agreement except for cytosine. Unfortunately, 
the experimental values were defined as the difference be- 
tween the vertical ionization potential and the adiabatic 
ionization potential, with values extracted from differ- 
ent sources [22|, [2(| |27|. As a result, in general the final 
values can be inaccurate due to the use of different ex- 
perimental techniques and agreement of the theoretical 
data with the experiment is not completely reliable. For 
the nucleobases, the large geometry relaxation between 
the neutral and oxidized states, particularly the change 
in C-C, C-N and C-0 bond lengths, is the cause of the 
large magnitude of the inner-sphere reorganization en- 
ergy. In the case of the pair formation, the inner-sphere 
reorganization energy is defined as the geometry relax- 
ation of both nucleobases of the pair and is decreased by 
the flexibility of the hydrogen bonds in the neutral and 
ionic geometries. The unit charge, spread between the 
two nucleobases in the pair instead of one, decreases the 
geometry relaxation of each nucleobase as well. For the 
A-T pair the inner-sphere reorganization energy is found 
to be one-third of the sum of the reorganization energy 
of the separated adenine and thymine. That is a result 
of insignificant geometry relaxation of the nucleobases 
itself during the oxidation process because of the high 
flexibility of the two hydrogen bonds (opening transla- 
tion [23j]). There are three hydrogen bonds between the 
nucleobases in a G-C pair, which restrict the translation 
and rotation flexibility of the nucleobases. This causes 
the not so significant decrease of the inner-sphere reor- 
ganization energy for the G-C pair compared to that of 
the A-T pair. 

Further, the hydrogen bonds are the channels for 
charge transfer between the nucleobases. In the oxidized 
state the hydrogen bonds participate in the charge trans- 
fer between the nucleobases to bring the pairs from the 
nonequilibrum state, where the charge is localized only on 
the nucleobase with a lower IP, to the equilibrium state, 
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1 2 3 4 5 6 



sequence number (n) *■ 

FIG. 5: The inner-sphere reorganization energy Ai and the 
vIP values versus the number of pairs in a DNA duplex 
oligomers (A-T) n and (G-C)„ performed with UB3P86/6- 
311+G* 



where the charge is spread over the base pair [32j]. We 
consider the A-T and G-C base pairs as a single state for 
the following calculations of \i . The stacking of the base 
pairs into the (G-C)„ and (A-T) n oligomers leads to a de- 
crease of the inner-sphere reorganization energy Aj and 
a decrease of the vIP as well. The results are presented 
in Figure [5j where the decrease of \ is seen to occur 
due to the contribution of the rotation and translation of 
the base pairs relative to each other and to the spreading 
of the charge between the pairs. According to our data, 
with elongation of the (A-T)„ and (G-C) n oligomers the 
twist of the base pairs mostly contributes to the decrease 
of the geometry relaxation of each nucleobase and in a 
reduction of the A^ . The decrease of the energies of the 
adiabatic IP (see Fig. [3]) and the inner-sphere reorga- 
nization energy (see Fig. [5]) provide the decrease of the 
vIP, which is the sum of above two components. 

As we mentioned above, A^ depends on the charge dis- 
tribution over the chain. The electrostatic potential dis- 
tribution in the (G-C)„ and (A-T)„ oligomers and respec- 
tively the residence of the HOMO in the oligomer centers 
provides the localization of the charge on the central gua- 
nines and adenines in the oxidized state. We have cal- 
culated the charge distribution as the difference between 
oxidized m and neutral states with Mul- 

liken population analysis (33|. In the (G-C)„ and (A-T)„ 
sequences the charge is distributed along the chain and 
is characterized by the low charge density at the DNA 
molecule sides. For example, the density of the atomic 
partial charge localized on the n=l site is lower than that 
at the chain center by 0.06 coul for the (A-T)„ and by 
0.25 coul for the (G-C)„ sequences. For the (G-C)4 se- 
quence our results are in agreement with the data in Ref. 

Therefore, the charge accumulation in the oligomer 
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sequence number (n) 



FIG. 6: The inner-sphere reorganization energy A;,d corre- 
sponding to the single base pairs within the (A-T) n and the 
(G-C)n oligomers, where n = 3 and n—5 are calculated with 
UB3P86/6-311+G*. 



centers in the oxidized state produces the maximum ge- 
ometry relaxation in the center of the DNA chain. Wc 
have performed an estimation of the geometry relaxation 
of the separated base pairs \™ D within the optimized 
geometries of the (A-T)„ and (G-C) n oligomers, where 
n = 1 ... 6. The simulation results of A™ D for n=3 and 
ri=5 are presented in Figure [5] Clearly, for the (G-C)„ 
sequences the difference of the structure relaxation at the 
sides of the chain and in the center is significant than that 
for the (A-T)„ sequences. The behavior of these curves 
repeats primarily the charge distribution in the (A-T)„ 
and (G-C)„ sequences. 

The difference between the inner-sphere reorganization 
energy of the (A-T) n and (G-C)„ oligomers should pro- 
vides the larger magnitude of the vibrational coupling 
constant for the G-C pairs than that for the A-T pairs, 
and larger for the guanine than that for the adenine (see 
Table IJ). 



IV. CONCLUSIONS 

We have performed accurate quantum-chemical calcu- 
lations to determine the electron coupling and the inner- 
sphere reorganization energy for the (A-T)„ and (G-C)„ 
DNA oligomers, where n — 1 ... 6. The electronic cou- 
pling between the two neighbor nucleobases within the 
same strand decreases exponentially with increasing of 
the base pairs number n participating in the chain for- 
mation. The 7i >4 is the sequence number required for 
an accurate evaluation of the electron coupling in the 
DNA molecule. The orbital distribution in oligomers 
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with the HOBO residing on the central nucleobase have 
been found to be the main reason for charge accumulation 
on the base pair located close to the chain center. The 
charge distribution in the chain determines degree of the 
the geometry relaxation of the base pair during the oxi- 
dation process in dependence on their location within the 
oligomer. Therefore, the base pairs in the chain center 
have stronger geometry distortion during the oxidation 
process. Such results are in good agreement with the 
theory of polaron formation in the DNA molecule, where 
the maximum structure distortion occurs in the polaron 



center (35j . 
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